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ABSTRACT 

The Electrochemical Machining ( ECM) is one of the most 
widely used unconventional machining process. However, the 
potential of this. highly capable process has not been fully 
exploited due to the inherent complexities of the process, 
like , simultaneous occurrence of two phase fluid dynamics, 
unsteady state heat transfer, mass transfer, thermodynamics, 
electrostatics and electrochemistry between moving boundarieSo 

1 

This complicates the problem of anode shape prediction in 
ECM many folds. Majority of the models available for anode 
shape prediction are based on simplified assumptions. 

In the present work, a modified anode shape prediction 
model has been proposed which predicts the anode profile in 
stagnation, front, transition, side and stray current zones. 

This model accounts for the variation in different process 
parameters during machining. Experimental data has been used to 
ascertain the accuracy of the model. Comparison of analytically 
predicted and experimentally obtained anode profile reveals a 
reasonable agreement between the two. The effects of the 
assumptions made and the experimental errors and their possible 
contribution in causing the deviations between theory and 
experiments have been discussed* 
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Conventional ECM tool design is not accuiate enough to 
account for diff erent . complex phenomenon taking place in the 
domain of interest, A method has been suggested for determi- 
nation of the shape of the tool (cathode) (specifically 
accounting for the effects of stray current zone) to produce 
the required work ( anode) shape, while machining under 
specified conditions , The tool design model has been further 
extended to predict the tool shape which will produce a spike 
of specified size in the stagnation zone. This model can be 
used for analysing the problems of tool design for external 
shaping by ECM, This method accounts for variation in different 
process parameters. The cathode shape prediction problem has 
been attempted for electrochemical drilling with a bare tool. 
Tools have been designed for producing profiles obtained 
experimentally. Tool design has also been carried out for 
straight sided holes with the top surface flat, as well as for 
a reasonably small spike .in the stagnation zone. Comparison 
of designed and experimental tool shapes reveal good agreement. 



1 


CHAPTER I 

INTRODUCTION AND LITERATURE SURVEY 

1.1 ELECTROCHEMICAL MACHININ3 ( ECM) 

Advanced technological developments in industries like 
aerospace, nuclear etc, have been accompanied by development 
of materials which are exceedingly difficult to machine. These 
materials have been developed to meet the demand of high strength 
and heat resistance. These rapid developments in the field of 
materials have given an impetus to the manufacturing technology 
to develop, modify and discover newer technological processes 
with a view to achieve results that are far beyond the scope of 
existing conventional manufacturing processes. This has led to 
the development of new techniques of machining such as ECM, EDM, 
USM, AJM, EBM etc. These unconventional processes are capable of 
providing effective solutions to the problems imposed by the 
increasing demand for accurate machining of high~strength- 
temperature-resistant ( HSTR) alloys, the requirements of parts 
with complex geometries, and materials so hard as to defy machin- 
ing by conventional methods. 

Of all the unconventional machining processes. 
Electrochemical machining is the one that is most widely used, , 
ECM is a process of metal removal from an electrically conducting 
workpiece by controlled and accelerated anodic dissolution. 

High velocity electrolyte flows between anode and cathode 
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(Figure 1,1) subjected to a small d,c, voltage and a small 
interelectrode gap (IBG).is maintained between them. The 
flowing electrolyte helps not only in allowing the high rate 
of metal dissolution, but also carries away the reaction 
products and heat ge*nerated during the electrochemical 
reaction. 

The advantages of ECM process lie in its ability to 
machine very hard metals without causing any tool wear, 
high metal removal rate, no structural damage to. the workpiece 
and capability to machine complicated shapes with good surface 
finish. The various areas of application [ 1 1 of ECM are drilling, 
boring, turning, cavity sinking, milling, grinding etc. Although, 
the applications of ECM are many, its capabilities have not 
been fully exploited due to lack of clear understanding of 
complex nature of the process of metal removal and non-existence 
of an efficient tool design methodology,. 

The metal removal rate (MRR) in ECM depends upon a large 
number of interrelated parameters. Most of the machining 
parameters such as electric field distribution, electrolyte 
condbctivity ( K) , temperature rise (At), valency of dissolution 
(Z) , electrolyte flow velocity (U), etc, are varying simultane- 

ously and continuously along the electrolyte flow path with 
machining time. Further, the interaction of these parameters is 
of complex nature which makes it further difficult to understand 
the exact mechanism of metal removal. Moreover, the electrolyte 
in the IB3 is a mixture of liquid electrolyte, gases evolved 

{ 
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during electrochemical reaction and precipitates. In majority 
of cases, it is difficult to evaluate the flow parameters due 
to the effects of electrolyte starvation and cavitation. Thus, 
ECM is a fairly complicated phenomenon since it involves 
simultaneous occurrence of two-phase fluid dynamics, heat 
transfer, mass transfer, thermodynamics and electrochemistry 
betv/een moving boundaries. 

Further, the accuracy of machining is governed by how 
accurately the tools can be designed, wbich in turn depends upon 
the accuracy of anode shape prediction model. With this in 
view different models available in literature for cathode and 
anode shape prediction have been discussed below. 

1 .2 TOOL DESIGN FOR ECM 

The tooling problem is of great concern because it is 
usually linked with high cost, A direct method for deciding 
upon the tool shape is almost impossible for all but the 
simplest of work shapes. In practice, a 'trial and error* 
method usually has to be relied upon, with repetitive machin- 
ing runs having to be made to obtain the correct tool geometry. 

The designing of tools for ECM can be classified under 
two heads. First is the determination of shape of tool 
together with optimum machining conditions necessary to 
produce the required work shape [ 2 ] > The second aspect 
is a practical one, which is concerned with making the tool 
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of an appropriate material, fixing it in the machine, 
connecting it to the power supply and passing an adequate 
supply of electrolyte between tool and workpiece. The 
design requirements for some of these aspects might turn 
out to be conflicting and therefore may call for certain 
modifications in the tool geometry. Here, only part of 
the first aspect of tool design has been considered, i.e,, 
designing the tool shape for producing desired work profile, 
ECM deals with computation of tool shape which under specified 
machining conditions would produce a work having a prescribed 
shape or profile. Inversely, the problem of anode shape 
prediction deals with determining the work shape obtainable 
from a tool of known geometry when machining is performed 
under specified machining conditions [4,5], The workpiece 
profile is usually computed from the equilibrium gap. The , 
lEG is not uniform in ECM and it is influenced by a number 
of parameters such as electrolyte flow rate, work material 
microstructure, change in valency of dissolution of the 
work-material during machining, type of electrolyte, current 
density, stray current attack etc. The influence of all 
these continuously varying parameters is not well known, 
which makes the anode shape prediction difficult. 

In general, the practical methods like ’trial and 
error’ methods for obtaining the tool shape are expensive, 
time consuming and inaccurate. Several attempts have been 
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made to develop a comprehensive model of the mechanics of 
metal removal in ECM, Some of the models for anode shape 
prediction and cathode shape design are discussed below. 


1 .3 COSO METHOD 

This method was first proposed by Tipton [6,7] for 
the computation of equilibrium gap for any given operating 
condi-tions. This method excludes the considerations of the 
mode of electrolyte flow, overpotential, variation in electrolyte 
concbc tivity, heat transferred to the electrodes etc. This 
method can be applied to the electrode regions where the 
electric field can be assumed to be normal to the electrode 
surfaces. In this method equilibrium workshape is computed 
corresponding to the tool whose profile has to be approximated 
by several planar sections inclined at different angles, say, 

Q as shown in Figure ( 1 ,2) , 


For plane parallel electrodes with their surfaces 
normal to the direction of feed, the equilibrium gap, Y^, is 
given by 


= 


K (V V) 


F P 


w 


( 1 . 1 ) 


where, K is the electrolyte conductivity, (V- V) is the 
effective voltage, F is the Faraday's constant, is the 

mass density of the work and f is the feed rate of the 


tool 
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The equilibrium gap between any section of cathode 
surface and the corresponding surface of the anode which 
is parallel to it is given by, Y^/Cos0o The angle '0' is 
measured between normal to the anode surface and the 
direction of the cathode feed. The cathode shape for a 
given anode surface can be computed from the following 
equations , 


Let Y = f ( x) represent the anode surface. Then any 
point, sayA(x, y) , on the anode surface has a corresponding 
point, say B(x^»y^) on the cathode surface. The gap width 
between A and B then becomes Y^/CosO (Figure 1,2), We have, 

Y-Y^ = AB Cos© = Yg ( 1 ,3) 

and 


X.| -X = AB SinO = Y^ tan©, 
or 

By substituting the values of X and Y from equations , 
(1.3) and (1,4) in equation (1,2) gives the equation for 
the cathode surface. 


The ’Cos©’ method is generally applicable to the 
regions where the local radii of curvature of the anode and 
cathode surfaces are large compared with the equilibrium 
gap. This method excludes the effect of different process 
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parameters on predicted anode profile. Workpieces having ' 
sharp corners and complex shapes cannot be analysed with 
this method, 

1 .4 N0M03RAPHIC APPROACH 

K'dnig [4] and Ganesh [s] suggested the nomographic 
approach to evaluate the equilibrium anode shape. K'dnig [4] 
has prepared a nomogram for calculating the side gap for 
known values of equilibrium gap, tool radius and bare tool 
length. Ganesh [Sjhas developed separate nomograms for 
evaluating the equilibrium gap and tool feed rate. Tool 
feed rate is calculated for the given current and area of 
cross-section of the tool. The equilibri.um gap is then 
evaluated from the thus calculated feed rate and the given 
voltage. Nomograms facilitate planning of ECM operations, 
but since they have been prepared with many simplified 
assumptions they have not gained much popularity. Moreover, 
they are a pplica ble for specific machining conditions and 
work-tool combinations,, 

1 .5 CONTINUITY METHOD 

Lawrence [ 9 ] developed a graphical technique for 
tool design by solving the field equation in the gap. 

But this procedure gives low accuracy and is tedious. He 
subsequently adopted numerical methods to solve the field 
equation using digital computer. In this method, starting from 
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a defined equipotential work shape boundary, a model of 
electric field is constructed, where each new equipotential 
surface is a possible tool shape which will produce the 
original work shape. This method rests on the assumption 
that Ohm's law applies across a small element of electrolyte 
and that the current flowing between a pair of flux lines in 
an electric field is constant. It does not take into account 
the variation in electrolyte properties in the IB3 along 
the flow path, 

1,6 CATHODE SHAPE PREDICTION BY DRAWING CURVILINEAR SQUARES 
IN THE IBS 

Tipton [10 ] has developed a method for determining 
the tool shape by drawing curvilinear squares between the 
work and tool boundaries. The work surface is defined numeri- 
cally in XY plane, Y is taken as feed direction and coordinates 
are spaced by constant X increment equal to a fraction of 
the equilibrium gap. The work boundary is considered as 
equipotential with each coordinate representing the end of 
a current flow line. In equilibrium, the normal current 
density being constant, a constant current increment Aj 
flows between each of the adjacent coordinates on the work 
boundary since they are equally spaced in X-direction. By 
using each adjacent pair of coordinates, it is possible to 
draw the curvilinear squares and to construct the first row ' 
of potential field distribution between work and tool as shown 
in Figure 1.3. If the X coordinate spacing is Y^/n, then 
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applying the above process repetitively will result in 
th 

n equipot ential which is also the tool surface. Each 

intermediate surface is also a possible tool surface for 
equilibrium gaps as Y /n, 2Y /n, 3Y /n, 4Y /n, etc. 

This method is error prone, especially wAien the 
process is repeated several times, since it is based on 
differentiating the input data for drawing curvilinear 
squares, and when repeated any random error in the input 
data gets magnified. Complete failure of calculation 
occurs when two current flow lines cross each other. In 
such a case, the previous equipotential could be used as a 
tool with the fractional equilibrium gap associated with 
it. This model does not account for variations in the 
process parameters, 

1 .7 NUMERICAL METHODS 

The emergence of fast digital computers has led to 
the application of numerical techniques to simulate the ECM 
process. In ECM, during machining, local current density 
causes the change in the shape of workpiece. Further, 
the current density in turn changes with the changing work- 
shape, To determine the current density, it is necessary 
to solve Laplace equation for the field distribution in 
the IBG. The three popular methods used to determipe the 
potential distribution are discussed beloWo 
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1 .7.1 FINITE DIFFERENCE TECHNIQUE (FDT) 

Tipton [6 land, Hopenfeld and Cole [11] applied 
FDT to determine potential distribution in the IBS, 
assuming constant conductivity and temperature. 


Figure (1.4) shows the tool work surfaces drawn 
in square meshes containing a set of grid points of general 
coordinates ( i, j). The grid points on the cathode boundary 
are assigned zero potential values, and those on the anode 
boundary are set to some known values. The initial potential 
at the grid points within the IB3 region are set by linear 
interpolation along the vertical grid lines between anode and 
cathode. These values must be progressively adjusted until 
they satisfy a finite difference equation corresponding to 
Laplace's equation in the region between the boundaries. 

For a point 'O' located in a mesh of spacing h(i,j), 
the finite difference equation corresponding to Laplace's 
equation becomes. 


= 

3 X 3 y 


(lo5) 


Using equation (1.5), the potential at 'O' can be 
obtained as. 


^^i+1.J ^i-1.J ^i.i+1 ^i.i-1? 


^ i, j ) 


4 


( 1 . 6 ) 
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In some cases, all points on the anode and cathode 
boundary ma y not lie on the grid points. In such cases, 
it is customary to relax the potential by linear inter- 
polation, Nanayakkara and Larsson [l2] have suggested 
the use of irregular grids along with regular grids to 
account for complex geometries, 

1,7.2 FINITE ELEMENT METHOD (FEM) 

In majority of cases, the IB3 is complex shaped, ^ 

As a result, the cathode and anode boundaries cannot be 
matched accurately using square meshes, which introduces 
further approximations. In order to overcome this problem, 

Jain and Pandey [13] suggested the application of FEM to 
analyse the complex shaped lEG, 

Jain [3] has developed FEM models, FET-11 and FET-22, 
for the prediction of anode profiles in case of plane parallel 
machining. The model FET-11 is based on use of Ohm's law 
for the calculation of current densities in the IB3, and FEM 
has been applied for one dimensional problem to obtain 
temperature distribution within the IB3, In the model FET-22, 
Laplace equation has been solved to obtain potential distribution 
within the IB3 from which current densities are calculated, 

Jain also extended model FET~11 (and named it as SGFET-11) for 
predicting the anode profile in case of rectangular and 
cylindrical deep hole drilling operations. All these models 
are capable to predict the anode shape accurately, by taking 
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into consideration, the variations in different process 
parameters. Murugan [l4] and Yogindra [15] have extended 
FEM models for anode shape prediction using bit type of 
tools. Ravi Raju [16] has modified model SGFET-22 for 
cathode design of bare and bare bit tools. In all these 
models however, the stray current zone and the stagnation 
zone have been neglected. 

1 .7.3 BOUNDARY ELEMENT METHOD ( BEM) 

Both the methods, FDM and FEM are domain based 
requiring frequent grid regeneration to compensate for 
the movement of work boundary during simulated ECM, Also, 
both FDM and FEM yield solutions corresponding to internal 
as well as boundary nodes, but in most of the cases, internal 
solutions are unnecessary when predicting the workpiece 
boundary and designing the cathode shape. In BBA unlike FDM 
and FEM, boundary a lone is discretized into segments. This 
considerably reduces the amount of input data and its prepara- 
tion time. However, if required the solution to internal 
points can also be obtained. 

■■ f 

Narayanan, Hinduja and Noble [l7] have applied BEM, for 
the anode shape prediction in ECM, Linear and quadratic 
elements were used to represent the boundaries. It was 
concluded that with far less computational effort, a 
comparatively better accuracy is obtained using BEM. However, 
for a clear understanding of the process solutions at internal 
nodes are also required. 
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1,8 ANODE SHAPE PREDICTION IN STRAY CURRENT ZONE AND 
STAGNATION ZONE 

1 .8,1 STRAY CURRENT ZONE 

The distribution of electric flux lines during ECD is 
such that some portion of the workpiece which is far away 
from the tool also undergoes dissolution, due to current 
lines acting over this zone (Figure 1,5), This effect is 
called as stray current effect and the zone over which the 
stray current acts is called as stray current zone. Stray 
current attack is undesirable not just because of poor 
dimensional accuracy but also because etched and pitted surfaces 
are obtained in this zone, as reported by Hoar and Mears [18], 
Rather than predicting the anode shape in stray current zone 
for a particular type of electrolyte, the trend in literature 
was observed to be of minimizing the stray current attack by 
choosing an appropriate electrolyte or by modifying the 
properties of a given electrolyte, Boden and Evans [19 ] 
showed that the effect of stray current attack could be 
reduced to a certain extent by bringing about passivation 
in the stray current zone by addition of small amount of 
another suitable salt to the electrolyte. They suggested 
that the additive could be carbonate in the case of machining 
of Iron and Nickle, which would form a sparingly soluble 
reaction product with the workpiece and hence bring about 
passivation by precipitation of insoluble products in the pits 
of stray current zone, NaClO^ was found suitable for machining 
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of cobalt. LaBoda and Mcmillan [20] has suggested the 
use of sodium chlorate as electrolyte for machining of 
any steel workpiece including high alloy steels because 
of its low ’throwing power* v\tiich ensures that the cutting 
is concentrated in areas close to the cathode, as is 
desired. 

Bannard [21 ] has pointed out that the extent of 
stray current could be controlled to a certain extent by 
use of a "thief -anode" (Figure 1,6) which steals away the 
current flux which is not- normal to the tool face. The 
disadvantage of this method is the difficulty in the 
manufacturing of the tool„ 

1 .8.2 STAGNATION ZONE 

Little literature appears to be available till date, 
discussing about the profile of the spike that is- form^ 
in ECD of blind holes. Larsson and Muzaffaruddin [22 ] 
assumed radial current density distribution to determine 
the spike hei^t. The purpose of their analysis was to show 
how the shaping performance of electrolytes is related to 
their electrochemical properties. They used the spike height 
as an index to demonstrate the effect of solution properties 
on its shaping performance rather than to predict the actual 
spike profile. The use of a simple tubular drill with an 
uninsulated bore constituted a very sensitive test for shaping 
performance, as the hei^t of the spike produced at the bottom 



15 


of a blind hole is shown to be dependent on both polarisation 
and current efficiency. They showed that the development of 
a current density dependent voltage at the electrode- 
electrolyte interface, and efficiency of an electrolyte are 
far from negligible in determining the shaping performance 
of an electrolyte. This is particularly so v\hen machining is 
in generating mode where there is a wide range of current 
density. They concluded that polarization plays an important 
role in reducing the shaping performance of an electrolyte and 
should not be ignored in any method used to predict work- 
shapes* Here, the assumptions made are far reaching, since 
the current distribution at the top of the spike will 
certainly not be radial. Also, the change in electrolyte 
conductivity due to heating effects and the development of 
gas bubbles in the electrolyte has not been taken into account, 

1 ,9 PRESENT WORK 

From the above literature survey, it appears that no 
attempt has been made to determine the effect of stray current 
attack on predicted anode shape. Also, no analysis appears 
to be available to find out the profile of the spike that is 
formed in the stagnation zone djring electrochemical drilling 
of blind holes. For a comprehensive and meaningful analysis of 
anode shape prediction during ECD, the analysis of stray current 
zone and stagnation zone is very important. The work embodied 
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in the first part of this thesis has been aimed towards the 
development of anode shape prediction model to analyse the 
effect of stray current attack on the anode profile* The 
anode shape prediction model (93FEr-22) developed by 
Yogindra [15] has been modified by incorporating the stray 
current zone,' The model has been named as STZFET-22, This 
model has been extended to determine analytically, the 
profile of spike generated in stagnation zone during ECD, The 
analytical anode profile obtained by STZFET-22 has been compared 
with experimental profile obtained by Jain [3] « 

The second part of the work has been aimed towards the 
development of a conputer aided method for the tool design 
for ECM. The two dimensional cathode (tool) shape prediction 
model for ECD proposed by Ravi Raju [16] has been modified to 
account for stray current zone and stagnation zone. The anode 
shape prediction model STZFET-22 has been used in this design 
package for electrochemical drilling. Tools have been designed 
for producing the specified work shapes, while machining 
under prescribed machining conditions. Designed tool shapes 
have been compared with experimental tool shapes used by 
Jain [3] for electrochemical drilling. Design of electrolyte 
supply hole in the cathode has been ained to minimize the spike 
dimensions, especially the height. This is a case similar to 
the tool design for external shaping by ECM process. The tool 
design model has been named as STZDES-22ii 
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FIG, 1.1 BASIC SCHEME OF ELECTRO-CHEMICAL MACHI- 
NING PROCESS. 




(a) Principle of the cos6 method (b) Cathode design by cos 0 methods 


FIG. 1.2 Cose METHOD 





Fig. 1.3 The first step in calculating a tool shape by cal- 
culating the field distribution between work and 
tool as a series of curvilinear squares. 



X 


a) lEG-tool work system discretized (b) A squre mesh of spacing 
in square meshes h with general point 0(1, j) 


T 



S 


(c) Boundary forms irregular star with one short arm QT 

Fig. 1.4 Finite difference technique. 










CHAPTER II 


THEORETICAL ANALYSIS 


2.1 INTRODUCTION 

The ability to predict 'the variations in process 
parameters for any given operating conditions is a pre- 
requisite for proper design of tools in ECM, In this 
chapter, the expressions for various process parameters 
such as electrolyte flow velocity, electrolyte conductivity 
current density, temperature distribution etc, have been 
derived for electrochemical drilling. In the present work, 
ECD with outward mode of electrolyte flow has been analysed. 
However, the same analysis can be applied for the case of 
inward mode of electrolyte flow with small changes in the 
computational scheme. To make a comprehensive analysis, all 
the five zones (Figure 2,1) namely, stagnation, front, transi 
tion, side and stray current have been analysed individually 
with suitable assumptions mentioned therein. The stagnation 
zone has been divided into three zones as shown in Figure 
(2,1). The three zones are named as stagnation transition 
(S.T, zone), stagnation side zone (S.S, zone) and stagnation 
stray current zone (S.S.C, zone). 
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2.2 ANALYSIS OF ELECTROCHEMICAL DRILL INS 
2.2.1 ELECTROLYTE FLOW VELOCITY (U) 

In the present work, the mode of electrolyte flow 
is outward and separate expressions have to be derived for 
stagnation, front, transition, side and stray current zones. 
Th e a ssumptions made for calculation of electrolyte flow 
velocity are: 

C i) Electrolyte is an ideal fluid. 

( ii) The volume of electrolyte flowing at any section 
per unit time is constant. 

( iii) Cavitation and starvation of the electrolyte are 
neglected. 

( iv) El ec tnelyte velocity across the flow direction is 
assumed to be uniform 


2.2.1 .1 STAGNATION ZONE 

As the inner tool comer radius is very small 
electrolyte flow in S.T. zone, S.S. zone and S.S.C. zone 
(Figure 2.1) can be considered to be in axial direction. 

The following analysis to determine electrolyte flow velocity 
holds good in all the three zones, i.e. S.T, zone, S.S. zone 
and S.S.C. zone. Referring to Figure (2.2), 

Ar = (i^r^ - -^(r^ - Y)^) 


or. 


= %i2T^ Y - Y^) 


( 2 . 1 ) 


and 
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U 


n 


( 


Q 

2T 

y ~ y ) 


where, Q is discharge of electrolyte/unit timso 


( 2 . 2 ) 


t 


2.2.1 .2 FRONT ZONE 

In this zone, the electrolyte flow is assumed to be 
in radial direction. Electrolyte flow velocity at any 
radius (Figure 2«3a) frcxn the axis of the tool is given by, 



Q ^ 
2 'n: r V 


(2.3) 


2,2.1 .3 TRANSITION AND SIDE ZONES 

For calculation of electrolyte flow velocity in 
transition zone, following assumption has been made. If 
the X coordinate of the node at which flow velocity is being 
calculated is less than the outer radius of the tool ( r 2 ) 
then that particular node is considered to be in front zone, 
otherwise it is considered to be in side zone. In side zone, 
the electrolyte flow is assumed to be in axial directi on(i: 


From Figure (2,3c), 


= 71 ( r2 + Y)^ - r2 

or 

= ii(Y^ + 2r2 Y) 


and 
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U = « 2 (2o4) 

7C ( + 2x^) 

2.2.1 ,4 STRAY CURRENT ZONE 

For finding out the velocity of electrolyte, Laplace 
equation has been solved for velocity potential. Here, 
the following assumption has been made. In stray current 
zone at edge AD (Figure 2.4), velocity is perpendicular to AD, 
which means vj5 = 0, i.e, = constant along edge AD, For 

convenience, jS has been taken to be equal to zero along edge 
AD. 

t 

For a two dimensional case the Laplace equation can 

+ = 0 . 

9 y 

Using simple triangular elements (Figure 2.4), 
the velocity potential (0) has been assumed to vary linearly 
within the element throughout the. solution domain. This 
velocity potential (0) can be represented by a polynomial 

0® = a + bx + cy , ( 2.5) 

The equation ^,5) can also be written in terms of the 
interpolating functions, N^, at each node, 


be wr itt en as, 
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i.e. jO® = + n^S (>2 + ^ 2^3 = 10 

Substituting equation (2,6) in equation (2,5), 
the residue becomes 


( 2 , 6 ) 


R® = 








a y 


By Galerkin’s method. 


// dx dy = 0, i = 1,2,3 

i,e. 

// N.. (-^ + — 1“) dx dy = 0 , 

3 X 3 y 

Integrating by parts w,r.t. x for the first term and 
w,r,t, y for the second term, gives 

( _J: + — i — ^ ) dx dy = 

X X y Y ^ 

(If" n, - |f V 

Where, n and n are the direction cosines along X and Y 
axis respectively, on the boundary of the domain. 

Using equation (2,6) and writing the equation ( 2,7) in 


matrix f orm gives 
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3N. 
fl (-i 
^3x 


3N I 
3 x J 


* !!!i pN 




dx dy 


N. V 
*^18 n 


\Miere, is the normal velocity. 


i.e. 


{0> 


= iv 

n 


( 2 . 8 ) 


where. 


K 


(e) 




and 




(e) _ 


^iB ^n 


dB, 


For an element, the matrix form of finite element 
equation becomes. 


wh ere. 


^11 

Ki2 

3 

K 21 

^22 

K23 

K31 

K32 

^33 

^ij 

= // 

3N^ 
^ Vx 


- 

02 

II 

1 

"1 


^3 


CO 

c 

> 






BN. 3N. 

i J. 

3 y 3 y 


)dx dy 


( 2 . 9 ) 


( 2 . 10 ) 


dB 


and 
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The interpolating functions for an element are 


^ (X2Y3 - X3Y2) + XlYj-Yg) + V(X3 - Xj) 

and 

N3 = 3 + ’^(^2 - >' l ) 

2 A 

where, 

A® = 5 (^1^2 - ^ <V3-W + • 

( 2.11 ) 

( q) 

After finding out element stiffness matrices [ Kj^ 
and normal velocity vector for each element, they are assembled 
for the entire domain to get a global stiffness matrix and a 
global velocity vector respectively. On applying boundary 
conditions, given by equation (2,12), 

0 along solid surfaces 

0 along open edge (2,12) 

= Velocity in the gap at 
the end of side zone 
(Figure 2,4) 


_ 

an - 
0 = 


a n 


the system will yield a set of linear equations, which could 
hp s (-lived for ^ bv anv standard techniaue. Here, Gauss 
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elimination technique has been used. 

Once the velocity potential (0) is obtained at each 
node, the velocities in stray current zone are calculated 
using equation (2,13). 

U = u^ + v^ (2.13) 

where, 

2 A ® 

0l(X2 - X3) + 02(X - X^) + 03 (X^ - X^) 

2A 

0 

and A is given by equation (2,11). 

2.2,2 CURRENT DENSITY (J) AND CONDUCTIVITY (K) 

The current density can be determined either from 
Ohm's law or from potential distribution obtained by solving 
Laplace equation. 

The application of a particular case depends upon 
the nature of interelectrode gap (IBS) and electric potential 
lines. Ohm’s law can be applied only when the current lines 
are straight and normal to the electrode surface. But this is 
not true when electrodes with small radii or complex shapes 
are used. In such cases, the current density calculated from 
potential distribution obtained by solving Laplce's equation 
gives more accurate results. The current density at any point 
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on the workpiece is calculated as 

J '= K ^ (2J4) 

The electrolyte conductivity, K, is a function of 
electrolyte temperature (T) as well as void fraction 
T.ipton [7 ] has proposed equation (2,15) for calculation 
of electrolyte conductivity, which takes into account 
the effect of electrolyte temperature and hydrogen gas 
liberation, 

K = K. (1 + a at) (1 - (2,15) 


where, 

is the initial electrolyte conductivity and a 
is the temperature coefficient of electrical conductivity 
of electrolyte. 

In equation (2,15), Hopenfeld and Cole [ll] assumed 
a value of void fraction exponent, n = 1 ,5 for the case of 
uniform void distribution. For the case of non-uniform 
void distribution, especially where bubbles are concentrated 
near cathode, Thorpe and Zerkle [23] suggested a value of 
n = 2,0, 

In stray current zone and stagnation zone, the effect 
of void fraction on conductivity is neglected, hence in 
these zones, the electrolyte conductivity is a function of 
temperature alone. 
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2.2.3 TEMPERATURE (T) 

Following assumptions have been made in deriving 
the equations for temperature distribution: 

C i) Total heat generated is only due to Ohmic heating. 

The heat generated by other sources such as friction, 
chemical reactions etc. is negligible. 

( ii) Specific heat and density of the electrolyte remains 
constant • 

(iii) Heat transfer through the electrodes is considered to 
be negligible, i.e. all the heat generated is retained 
by the electrolyte. 

(iv) The flow velocities are such that the frictional 
heating of fluid is negligible,. 

2 .2 . 3 . 1 ST AG NAT ION ZONE 

In S.T. zone, S.S. zone and S.S.C. zone (Figure 2.1), 
the electrolyte flow is assumed to be in axial direction. 

The following analysis for calculation of temperature is appli- 
cable to S.T. zone, S.S. zone and S.S.C. zone. Consider 
an element of length dx at a distance x from the starting of 
S.T. zone (Refer Figure (2.2)),, 

By Ohm’s law, current through this element is given 
by, 

I ^ K 2-Ti;(r^ - y) dx . (2.16) 
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Rate of heat generated in the element is given by 

^ . (V -AV)^ K 2Tc(r^ - y) dx 

dH = R. = y — • 

(2.17) 

Rate of heat carried away by the electrolyte, 

dH = ^2.18) 

Equating ( 2.1 8) and (2.17) yields 

(V -AV)^ K 27t(r^ - y)dx 
Q dT = y-- • 

Using K = ( 1 +aAT), the expression for temper- 

ature distribution beccmes, 

<ff - - A' dx (Z.I'S) 

1 -raCT-Tj) 

wh ere, 

, (V -h'/y 2jt( - y) 

“ y 4PeC, 

Integrating equation (2.19) with boundary condition, 
at X = x^gf, T = gives, 

X - Tj_ = ^ f[exp a A' (x - ^^^ef^^ - (2.20, 

2. 2. 3 .2 PROMT ZONE 

Consider an element of thickness ’ dr' at a distance 
r, from the axis of the tool (Figure 2.3d). 
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Using Ohm's law, the current I through this element , 
is given by. 


I = 


(V -AV) K 2Txr dr 
Y 


( 2.21 ) 


Rate of heat generated in the element is given by. 


dH = I Rg = 


( V K 2nr dr 

Y 


( 2 . 22 ) 


Rate of heat carried away by the electrolyte. 


dH = Q ^e •* 


(2.23) 


Equating (2.23) and (2.;22) yields 


Q ^ C dT = 
e e 


(V -AV) K 27tr dr 

— 


Using K = ( 1 + a . ^T), the expression for temperature 

distribution becomes. 


dT 

1 +arT - l7) 


= A r dr 


(2.24) 


where. 


2 71 : (V - AV )‘ 


The integration of equation ( 2;.24) with boundary 
conditions, at r = r T = gives. 
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T - T; 


- X? 

exp [A (—5 — )]- 1 ^ 


(2o25) 


2, 2, 3, 3 TRANSITION AND SIDE ZONES 

In these two zones, the electrolyte flow is assumed to 
be in axial direction. Consider an element of length dx at 
a distance X from the end of front zone, i,e, from 
(Figure 2.3e) , 


By Ohm's law, 

(V -AV) K 2n:(r2 + Y) dx 


I = 


y 


(2.26) 


Fran equation (2,26), the rate of heat generated in 
the element is given by 


o ^ (V - AV) K 2n( r^ + Y) dx 

dH = I^ = 

e 


T 


'(2,27) 


Equating equations (2,27) and (2,23) gives. 


Q dT = 

^ 0 e 


(V - AV)^ K 2'n:(r2 + Y)dx 

^ ^ — 


(2,28) 


or 


dT 


1 + aVT - Ti) 


It I 

= A dx 


(2,29) 


whe re , 
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«, 2TrK.(V - AV)^ (ro + Y) 

A = 1 — ^ (2o30) 

Q P e 

Integrating equation (2,^9) with boundary condition 
at X = ^ref» ^ ~ ^ref ^ives, 

T - T. = ^ tB exp [aA*’ (X - ] - 1 } (2.31 ) 

where, 

B = 1 +a(T^ef-^i) ^2.32) 


2, 2. 3. 4 STRAY CURRENT ZONE ' 

Consider an element of length at a distaixe r’ 
from the axis of the tool (Figure 2,5), 

By Ohm' s law, 

(V - AV)K 2 Tcr’ dr' 

I = j ( 2 . 33 ) 


From equation (2,33) the rate of heat generated in the 
element is given by, 


dH 


= I = 

e 


(V ~AV) K 2 r« dr' 


( 2 o 34 ) 



Rate of heat carried away by the electrolyte is given 


dH = QPq ^e (2,35) 
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Equating (2,35) and (2,34) gives, 

QP^C^dT = d^'. . . (2.35) 

After rearranging the terms. 


dT 


1 + a(T - T^) 


= r* dr* 


(2.37) 


where. 


IV 


2 K. (V -AV)‘ 


1 Q P e Ce 


(2.38) 


Integrating equation (2,37) with boundary condition, 
at r’ = R, T = T__ gives, 

o Z 


A ^^, _.2 „ 2 ^ 


- i £8., exp[a^ ( r* ^ - R^) ] - 1 } 


(2,39) 


where, 


Bi = 'I ^ "^i^ '• 


(2.40) 


In calculation of temperatures in Front, Transition, 
Side and stray current zones, initial temperature is taken 
to be equal to temperature at the end of stagnation zone^ 
say, T-, Similarly, initial conductivity of electrolyte is 
taken to be equal to the conductivity at the end of stagnation 
zone, say, K^, 
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2.2.4 VOID FRACTION (a^) 

As the electrolyte flows downstream within the machin- 
ing gap, due to electrolysis, hydrogen will evolve at cathode. 
This hydrogen evolution, usually reduces the effective conduct ivi 
of the electrolyte, so that the local anodic dissolution rate 
varies dovmstream until an equilibrium is achieved. 


Considering a control volume in IBG (Figure 2.6), using 
the principle of conservation of mass, Thorpe and Zerkle [23] 
arrived at the governing differential equation, and the 
expression for void fraction was obtained by solving it as 
f ollows: 


a 


V 


* » 

B X 

3 — 

1 + B X 


(2.41) 


wh ere. 



P 

e 


E. 


P g ^a 


t 


and 


* 

X 


2L 

LS 


for rectilinear flow. 


* X(r. + X/2) 

X = for radial flow 

LS r^ 


S = 


^ A 

-Tp-r- 




Here, E is the electrochemical equivalent and the suffixes < 
a,e and g stand for anode, electrolyte and gas respectively. 
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In the present case, the gas slip ratio o in the variable B 
is assumed to be unity, ioe,, velocity of gas is same as that 
of velocity of electrolyte* 

The effect of void fraction on conductivity has however 
been neglected in the stray current zone and stagnation zone 
due to the mathematical complexity of the problem, 

2*2o5 ELECTROLYTE PRESSURE (P) 

The electrolyte pressure [l] required to maintain 
the rate of electrolyte flow in a particular gap is given by 

P = Pl + P 2 (2*42) 

where, P-| is the pressure required to over come inertia force 
and P 2 is the pressure required to overcome viscous force* 

The expressions for P^ and P 2 are different for laminar 
and turbulent flows. The type of flow in the lEG can be 
determined from Reynold’s number (R^^)* given by, 

P u^ D 

R = -^—2 (2.43) 

n P e 

where, P g is the viscosity of the electrolyte and D is the 
hydraulic mean diameter which is given by, 

n =s 4 X Area of cross section. 

Wetted perimeter “ 

For a cylindrical tube, D = d, for a rectangular channel of 
depth Y, D = 2Y and for the case in which fluid is flowing 
through the annulus, D = ( d^ — ^ 2 ^ • 
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For the laminar flow ( < 2000) with a parabolic 
velocity distribution, the pressures and P 2 are given by 

2 

and 

for a narrow tube of diameter d. 

The flow is usually turbulent ( ^ 2000) in ECM, 
From Bernouli’ s equation. 


32 P ^ U X 
® e 


(2,44) 



t 


The value of P 2 is given by Blasius equation. 


0,3164 P ^ X 

D — ^ - . 

2D 


(2o45) 


The pressure P 2 » in turbulent flow is a function of 
Reynold's number. At high Reynold's number inertia force 
will be more than the viscous force which results in 
higher value of P^ » 


2,2.6 FEED RATE (f) 

In electrochemical drilling, as the machining progresses, 
the side zone becomes tapered and it no longer remains a zero 
feed rate case. For the calculation of lEG in the side zone, 
the feed rate vector normal to the work surface should be 
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considered, as given by equation ( 2o46 ) 

= foCose ( 2 . 46 ) 

where, 0 is the angle of inclination between the feed direction 
and normal to the work surf a ce . The values of ’0’ for each 
element in the side zone and transition zone are calculated, 
separatel y. 

If and are the interelect rode gaps in side 
zone at points 1 and 2 respectively ( Figure 2,7), then the 
angle 0 at the point 2 is given by equation (2,47), 

~1 ^ 2-^1 

« = tan ' (rr 7^) 4 (2.47) 

^2 “ 

In the transition zone, the angle © at any point is 
given by the intersection angle between the normal drawn at 
a point and the tool feed direction (Figure 2,7), 

In stagnation zone analysis, equations (2,46) and (2,47) 
have been used for calculating the feed rate in stagnation-side 
zone (Figure 2,1), Since stagnation- transit ion zone is similar 
to transition zone, feed rate in stagnation-transition zone has 
been calculated on the same lines as feed rate in transition 
zone. Unlike the other zones, in stray current zone feed rate 
does not play a role in determining the anode profile and is 
hence ignored. 
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2.2.7 IMTERELECTRODE GAP (lEG) 

Electrochemical, reactions would attain an equilibrium 
only after 't' approaches infinity and IB3 can be computed 
from equations (2,48) and (2,49) while machining with zero 
feed rate and finite feed rate [ 1 ] respectively.;, 

Y = (Yq + 2Ct)^/^ (2.48) 

where, 

C 

and 

Y - Y 

t = ^ [(Y^ - Y) + in (y^rT^ ) ^ ‘ (2.49) 

Equation (2,49) is an implicit equation and can be solved 
by iteration process. Equation (2;;48) predicts an infinite 
value of the gap as t approaches infinity. However, in 
practice, as the lEG increases, the current decreases, so 
that the MRR gradually decreases and for a very large gap 
the process would come to a stand still.,; 

For a very small interval of time At, Jain and Pandey [15 
derived an equation (2,50) for determining lEG which holds good 
for both, zero and finite feed rate , 

Y = Yq + (C« - f ) A t (2,50) 

where. 
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r.’ TiJ A 

^ p FZ ■■ 

^ w 

Equation (2o50)'is based on the assumption that the 
current density ’J’ remains constant over a small element of 
length, dX and for a small interval of time, At. This 
equation yields values of Y which are very close to those 
obtained from equations ( 2o48) and (2,49) and agree well with 
experimental data [3] , In the present work equation (2,50) 
has been used. 


2.2,8 CURRENT EFFICIENCY (h) 

In deriving the expression for metal removal rate by 
Faraday's laws of electrolysis, current efficiency was assumed 
to be 100?4, But in actual machining process, the MRR depends 
on the current efficiency achieved [ 1 ] , i,e,, the proportion 
of total current that is used for the removal of metal from , 
the anode. The current efficiency is given by. 


Actual metal removed 
Theoreticai metal 
removed 


X 100% 


(2,51) 


In actual practice, the values of current efficiency 
may be higher or lower than 100% depending on the conditions 
of machining. Current efficiency values less than 100% are 
obtained because of the side reactions at the anode, e.g,, 
evolution of oxygen gas or Fe ions instead of Fe ions 


etc 
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If the dissolution of the anode occurs at a valency 
lov/er than the one used for calculation, the current ' 

efficiency obtained may be higher than 100?^o Moreover, 
selective dissolution may occur at the material irregulari- 
ties such as at grain boundaries [24] and hence may cause 
removal of chunks of material© In this case, the actual 
material removed is more than the theoretical one, which 
results in current efficiency higher than 100%, 
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Tr. Zone — Transition Zone 

S.T. Zone —Stagnation - Transition Zone 

5.5. C. Zone-Stagnation - Stray Current Zone 

5.5. Zone— Stagnation - Side Zone 


Fig. 2.1 Various zones in ECD 



(a) Electrolyte Flow Area (b) Current Flow Area 

Fig. 2.2 Electrolyte and current flow area in stagnation 
zone tor outward mode of electrolyte flow. 
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(b) 



Froni zone (section at A-A) (c) Transition and side zone 

Electrolyte Flow Area 



— 2Tr(r2 +Y) — 



(d) Front Zone (e) Transition and side zone 

Current Flow Area 

Fig. 2.3 Electrolyte and current flow area in front, 
transition and side zones. 






Boundary Condition 
AD 0 = 0 

AB, CD (ad/3n)=0 

BC (3<^/8n)=Vn 


Fig. 2.4 Finite element discretization & velocity boundary 
conditions in stray current zone. 



Fig. 2.5 Current flow area in stray current zone for out- 
ward mode of electrolyte flow . 






CATHODE 



Fig. 2.6 Variation in gap width along electrolyte 
flow direction [23] 



Fig. 2.7 Feed rate components in side and transi- 
tion zones [153 
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CHAPTER III 

■ tool design 

3.1 INTRODUCTION 

Electrochemical drilling is a process of making holes 
of any shape and size using the principle of controlled 
anodic dissolution. The electrolytic cell for an ECD process 
is illustrated in Figure (3.1). The cathode is simply a 
straight tube, shaped for obtaining the desired work geometry. 
The electrolyte is carried to the drilled hole in the work- 
piece (anode) through the tool, and metal removal takes place 
as the electrode is fed in the downward direction. 

3.2 PREDICTION OF ANODE SHAPE 

For the purpose of analysis of the process, the entire 
machining region has been divided into five zones (Figure 3.2), 
namely, stagnation, front, transition, side and stray current 
zone. To reduce the complexity in the present problem, the 
stagnation zone has been analysed separately from the other 
four zones. Further, depending upon the location of inlet and 
outlet of electrolyte, the electrolyte flow mode can be classi- 
fied as outward, inward or side flow. 

In the present work, the finite element technique has 
been applied to the circular hole drilling with bare tool for 

, T 

outward mode of electrolyte flow. Anode shape prediction 
model developed by Jain and Yogindra [15] has been modified to 



48 


predict the anode shape in all the five zones mentioned 

above. In the anode shape prediction model developed 

here ( STZFET-22) , the field distribution is obtained by 

solving the Laplace equation = 0) with 

3 X ay 

the following boundary conditions: 

t 

(i) jQ = 0 along cathode (tool) 

( ii) 0 = (V - AV) along anode (work) 

(iii) = 0 (along open edges),. 

Since ECD is a moving boundary problem, the domain 
shape and size will change after each computational machining 
cycle time. At sec. As the tool penetrates in the work surface, 
the total number of nodes would also increase. The coordinates 
of each node, therefore should be recomputed after each compu- 
tational cycle. Figure (3,3) shows the finite element discre- 
tization and renumbering of nodes after one computational cycle 
in all the zones except stagnation zone. Figure ( 3 , 4 ) 
illustrates the same for stagnation zone analysis. Here, an 
assumption has been made that the feed to the tool is given 
only after machining has taken place for a time per5od of 4 t 
seconds, i.e, the process is considered as intermittent cutting 
rather than continuous one, as in actual ECM. 

Due to complex nature of ECM process, majority of 
models proposed for anode shape prediction in ECD, do not 
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account for the simultaneous variations in the process 
parameters such as temperature, conductivity etc. It is 
felt that for an accurate prediction of anode profile, 
the variations in all these process parameters with machining 
time and their effect on anode profile should be considered. 
These aspects have been included in the present model, 
STZFET-22. Some problems have been attempted using STZFET-22 
and the predicted anode profiles have been compared with 
experimental anode profiles, 

3.3 TOOL (CATHODE) DESIGN FOR ECD 

The problem of designing a tool to produce a required 
work shape is probably the single most important problem in ECM 
and in practise, methods of design are still empirical even 
to the point of ’trial and error'. The lEG between tool and 
work is a function of various machining parameters. Hence, 
the tool shape is not congruent with the work shape in 
the strict sense. The ECM tool design depends essentially 
on the accuracy of the corresponding anode shape prediction ^ 
model. Therefore, for accurate tool-design is ECM, the anode 
shape prediction model must account for the variations in 
process parameters while arriving at an anode shape. 

Reddy [25] has suggested a ’Correction Factor Method' 
for tool design in ECM. In this chapter, a two dimensional 
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model for predicting a cathode configuration which can 
produce specified anode shape under prescribed machining 
conditions has been described. This model has been sugges- 
ted by Ravi Raju [16 j and has been modified here to include 
stray current zone and stagnation zone. The anode shape 
prediction model (STZFEr-22) described earlier has been used 
for tool design in ECD, This model has been suitably 
modified for the prediction of work-shape with a tapered tool, 

3,4 GENERAL DESIGN PROCEDURE 

Tool design procedure for producing the required work 
shape while machining under given conditions is summarized 
below: 

STEP 1 - INITIAL TOOL SHAPE 

The first step is to assume a tool shape depending on 

the required work shape. For majority of the cases, the tool 
shape can be assumed as congruent to the desired work shape, 

STEP 2 - WORK SHAPE PREDICTION 

The second step is to predict the work shape 
obtainable from the tool shape assumed in the first step. 

Out of the nany models available for anode shape prediction 
in ECM, the models based on numerical methods like FDM and 
FEM, can predict the work profiles reasonably accurately, 

STEP 3 - ERROR AND CORRECTION 

The third step is to compare the predicted work 

shape with the required work shape and to obtain deviations 



between them. If the deviation is more than the specified 
tolerance, correction required to minimise this error has to 
be calculated, 

STEP 4 - MODIFICATION 

The fourth step is to modify the present tool shape 
by applying the correction appropriately. Using this modified 
tool, again the anode shape is predicted. Steps 2,3 and 4 are 
repeated till the predicted work shape with the desired tolerance 
is obtained. 

This is an iterative process and the flow chart for 
this process is shown in Figure (3,6), Tool shape is refined 
progressively in each design cycle. Finally, it will give 

1 

the tool shape which produces the desired work shape. 

From the above mentioned procedure, it is evident that 
any inaccuracy associated with the anode shape prediction 
model will be reflected on the designed tool shape , 

3,5 TOOL DESIGN FOR ECD 

Here, an attempt has been made to design bare tool 
for electrochemical drilling. Here, a case of only cylindrical 
hole (straight sided and tapered both) drilling has been 
considered. However, the model can be applied with slight 
changes, to drilling of non-cylindrical holes as well. Using 
above mentioned procedure the tool design model developed by 
Ravi Raj u [1 6] has been modified to incorporate the tool design 
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for stray current zone and stagnation zone. Details of the 
above mentioned steps for the present model are as given 
below. 


3.5.1 INITIAL TOOL SHAPE 

Initial tool shape has been assumed as a cylindrical 
one, with its radius equal to 1 mm lesser than the required 
hole radius for the case of tool designed for uniform diameter 
hole drilling. In transition zone, tool corner radius is aghin 
assumed as 1 mm lesser than the required work surface radius 
in the transition zone. This 1 mm value has been chosen 
arbitrarily, and any other value could as well be chosen, 

3.5.2 PREDICTION OF ANODE SHAPE 

The anode shape prediction model discussed in the 
earlier section has been used for computing the anode profile. 
The IB3's at the beginning of a computational cycle get 
modified in the following cycle. Also, the tool occupies a new 
position after each computational machining cycle. The 
modified coordinates are assigned to the corresponding points 
on tool and workpiece, to obtain the correct work shape at 
the end of a computational cycles 

3.5.3 CALCULATION OF ERROR AND CORRECTION 

The coordinate system for the five zones namely 
stagnation, front, transition, side and stray current zones 
is shown in Figure (3,5). The system of coordinates has been 
so chosai that the same error calculation principle holds good 
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for all the five zones. The error and correction at each 
node are calculated as follows; 

ERROR (I) = OWS (I) - RWS (I) (3ol) 

CORR (I) = -ERROR (I) (3.2) 

where, 

"til 

CORR (I) = Correction to be made at I node. 

*bh 

OWS (I) = Obtained work shape at I node. 

The obtained work shape at any node in front, transition, 
side zones and, S.T, zone and S.S, zone of stagnation zone 
is represented by Y co-ordinate of that node. However, in 
stray current zone and S.S.C, zone of stagnation zone, 

the obtained work shape at any node is given by the X coordinate 
of that node*. 

RWS (I) = Required work shape at I node* 

For representing the required work shape at any node 
the same argument given above for obtained work shape in 
different zones holds goodo 

The correction need not always be equal to error in 
all design cycles* It can be calculated from the trend of 
the previous design cycle. It may sometimes consume less CPU 
time [ 1 6] . From the second design cycle onwards t the correction 
may be calculated as given below^; 
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CORRj^(I) 


CORRj^_^(I) X ERRORj^(I) 
ERRORj^^^(I) - ERRORj. (I) 


(3o3) 


where, suffix K stands for the present design cycle and 
suffix K-1 stands for the previous design cycleo 


However, for the sake of simplicity, in the present 
work, equation (3,2) has been used. 


3,5,4 MODIFICATION 

The tool that has been used in the previous design 
cycle is modified by applying correction to Y coordinates of 
the tool at different nodes. The new tool shape is given by 
equation ( 3o4) , 


NTS (I) = Y(I) + CORR(I) (3,4) 

wh ere , 

Nrs( I) = Y coordinate of new tool shape at the I 
node 

Y(I) = Y coordinate of tool in previous cycle. 

Correction at any node (CORR (I)) may be positive or 
negative depending on the deviation between the obtained work 
shape and required work shape at a particular node. 

Unlike plane parallel electrodes ECM, in ECD a node on 
the workpisce is exposed to different nodes on tool for different 
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times. Thus, if a tool shape at a node is changed, it 
will affect the work shape at many other nodes on its 
way. This in addition to the complex nature of ECM owing to 
electrochemical, thermal and hydrodynamic factors, makes the 
tool designing more difficult in ECD with bare tool. The 
former problem, to some extent, is eliminated in case of 
electrochemical bit drilling* 


( 



Elec 


rolyte 
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Fig. 3.1 Schematic diagram of BCD with outward mode 
of electrolyte flow. 



Fig. 3,2 A complete picture of various zones in ECD 
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FIG 



3.6 FLOW CHART FOR GENERAL DESIGN PROCEDURE. 
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CHAPTER IV 

RESULTS AND DISCUSSIONS 

To ascertain the validity of any analytical model, 
it is essent5al to compare the analytically obtained results 
with the experimental data. With this basis in mind, in 
this chapter, a compar5.son between the analytical and experi- 
mental anode profiles obtained during ECD with bare tools has 
been made. Also, the tools designed for producing the experi- 
mentally obtained anode profiles and the tools used during ^ 
experimentation have been compared. However, since experi- 
mental work profile in the stagnation zone is not available 
it could not be compared with the analytically obtained work 
profile in the stagnation zone. Further, because of this, 
part of the tool designed for obtaining a desired work shape 
in stagnation zone could not be compared with tool used during 
experimentation, 

Jain [3] conducted experiments using brass as material 
for tools and low alloy steel castings and low alloy steel 
forgings as nsterials for anodes (Table - 1). The machining 
conditions fbr ECD experiments are given in Table - 2, The 
shape and size of anode profiles were measured by preparing a 
cast replica of a drilled hole, and projecting it on a tool 
maker's microscope. The method by which anode corner radius 
was measured, has been described in detail in reference 
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For analytical computation of anode profile, modified 
equation (2,50) suggested by Jain [3] has been used, 

4.1 ANODE SHAPE PREDICTION 

Using the model STZFEr-22 discussed earlier, anode 
profiles for different machining conditions have been 
obtained, assuming constant valency of electrochemical 
dissolution of the anode and taking a mean value for the 
efficiency of machining (Table -3), Voltage fluctuation 
history (Table - 4) has also been used while computing 
the anode profiles. This would obviously predict more 
accurate anode profile unlike the case in which mean voltage 
is used, since during experimentation voltage was varying 
rather than reimining constant. The anode profile obtained 
for stagnation zone by model STZFET-22 has been superimposed 
on that obtained for other four zones so as to get the profile 
for all the five zones (Figure 3.2), These profiles for 
different machining conditions have been shown in Figure 

(4,1), 

It can be seen from Figure (4.1) that in majority of 
the cases, the electrochemically drilled holes have the shape 
of an elongated S(/).. The large overcut at and near the top 
surface of the workpiece is attributed partly to the stray 
current attack and partly to the fact that this section of the 
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workpiece has been subjected to electrochemical dissolution 
throughout the period of ECD, A good quantitative agree- 
ment between theoretically predicted and experimentally 
obtained anode profiles was observed (Figure 4.1) » The 
small deviations between the analytical and experimental 
anode profiles in front, transition, side and stray current 
zones are due to the following reasons: 

(i) The valency of dissolution plays an important role 
in determining the shape and size of the computed 
anode profiles. The present analytical model assumes 
some constant value of valency of dissolution. But in 
actual machining, the valency of dissolution may be 
greater than the assumed one [26] and hence the 
assumed lower valency would result in larger 
ccxnputed overcut. In the present work, valency of 
electrochemical dissolution has been assumed to be 2, 

But as reported in [26] for the conditions prevailing 
in present case of ECD, iron may dissolve with a 
valency of 3 instead of 2, This would reduce the 
overcut and thereby narrow the deviation between the 
experimental and analytical anode profiles, 

( ii) The electrolyte conductivity within the IB3 is a 

function of both temperature and void fraction as given 
by equation ( 2,15) , The value of exponent 'n', in 
equation (2,15), depends on the nature of distribution of 
voids in the fEG. For uniform void distribution. 
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Hopenfeld and Cole [11 ] suggested a value of n = 1 .5,^ 
and for non-uniform void distribution, especially 
when bubbles are concentrated near the cathode, 

Thorpe and Zerkle [23] suggested a value of n = 2,0* 

But in actual machining, it is extremely difficult 
to know the exact void distribution. Owing to this 
reason, an average value of n = 1 ,75 has been assumed 
for ccxnputing void fraction. This approximation 
might also have contributed to some error in calcula- 
tion of overcut. In stray current zone, the void 
fraction has been ignored, to keep the problem 
mathematically tractable, and hence the condicti vity 
is assumed to be a function of temperature alone. 

This may also lead to some var3.ation between analytical 
and experimental anode profiles in stray current zone, 

(iii) In some of the tests, due to sparking or some other 
unavoidable reasons, the machining had to be stopped 
and the IB3 readjusted. This would change both 

the current and voltage variation history and hence would 
lead to deviations between predicted and experimental 
anode profiles, 

(iv) In the present analytical model, no consideration 
has been given to the effects of grain size, grain 
boundary attack, selective dissolution of various 
alloy components, and other metallurgical aspects. 

This is expected to affect the value of computed 
overcut. 
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Since no experimental results are available about 
the spike profile in the stagnation zone, the analytically 
predicted anode profile in this particular zone could not 
be compared. However, the analytically obtained spike 
profile seems to be some what larger than the one that could 
have been normally expected. This discrepancy can be 
assigned to the following reasons: 

(a) Reason (i) mentioned above, 

( b) No consideration for the effect of void fraction on 
conductivity. Conductivity has been assumed to be a 
function of temperature alone, 

(c) Use of simple triangular elements in finite element 
discretization of the stagnation zone which may not 
fit exactly with the work or tool boundary, leading 
to a n error, 

4.1 ,1 COMPARISON BETWEEN THE MODELS STZFET-22 and SBFET-1 1 
From Figure (4,1), it is evident that the model 
STZFET-'22 gives better agreement with experimental anode 
profiles in front, transition and side zones, as compared 
to those obtained from SBFET-11, Also, the model SBFET-1 1 
is incapable of predicting the anode profile in the stray 
current zone and stagnation zone. The model SBFET-11 ns a 
one-dimensional model and is based on the assumption that 
the lines of electric flux are straight and normal to the 
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electrode surfaces. The current density in the IB3 is 
calculated from the Ohm's law. But in complex shaped 
workpiece, especially at curved sections, the electric 
flux lines are not straight (Figure 1,5) and hence 
Laplace equation has to be solved to get accurate potential 
distribution in the domain of interest. In the model 
STZFET-22, current densities at different nodes are calculated 
from potent-ial distribution obtained by solving Laplace's 
equation. That is why, anode profiles predicted by STZFET-22 
are more accurate than the profiles predicted by SBFET-11 , 

The accuracy of the results obtained by finite element 
technique also depends on the type of elements selected for 
discretization of the solution domain. The model STZFET-22 
uses two dimensional triangular elements, which fit better 
in the curved boundaries, as compared to one dimensional 
linear elements used in model SBFET-11, This also contributes 
to the better accuracy of model STZFET-22, However, use of 
isoparametric elements would further increase the accuracy 
of computed results, 

4,1 ,2 EFFECT OF TEMPERATURE AID VOID FRACTION ON CONDUCTIVITY 

Figure (4,53) shows the variation in temperature, and 
void fraction and conductivity along the electrolyte flow 
path distance. The magnitude of increase in temperature is 
not much. From Figure (4,5a) it can be seen that the effect 



67 


of void fraction on electrolyte conductivity is more signi- 
ficant than that of tanperature rise. The electrolyte 
conductivity decreases drastically with the electrolyte flow 
distance. For clarity, the variation in temperature and 
conductivity along the electrolyte flow path in stray 
current zone is shown separately in Figure (4o5b')o It 
should be noted that in stagnation zone and stray current 
zone, the effect of void fraction on conductivity is neglected. 
Hence, in these zones the electrolyte conductivity is a 
function of temperature alone, 

4.2 CATHODE DESIGN 

Figure (4,2) shows a comparison between tool shape 
designed by correction method and the tool shape used during 
experimentation for the case of ECD with bare tools, 
under different machining conditions. Here, the required and 
obtained work shapes are coinciding because the design procedure 
is continued till the deviation between the required and 
obtained work shapes lie within an accuracy of 10 p-m. There 
is a good agreement between the designed and experimental 
tool used under the same machining conditions. Kinks, just 
opposite to the unevenness in the experimental profile, are 
observed In the designed tool shape. This is particularly 
observed at the end of transition zone. 

Figure (4,2) also shows the tool designed by Ravi Raj u [16] 
using model DES-22, The tools designed by model STZDES-22 and 
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model DES-22 show a close agreement between the two,, 

However, the model DES-22 is Sncapable of predicting 
the tool shape by taking into consideration the stray 
current effect. 

Figure (4,3) shov;s the (Resigned bare tool (using 
model STZDES-22) for producing a straight sided hole with 
top surface of the workpiece square with the hole walls 
(1,6, the top surface of the workpiece is flat and has 
not undergone any damage due to stray current attack). 

The tool shapes designed by model DES-22 have also been 
plotted for comparison (Figure 4,3), Once again, the stray 
current zone has been neglected in model DES-22, 

Figure (4,4) shows the designed tool (w^th an uninsula- 
ted bore) for obtaining the desired spike profile. The desired 
spike profiles for different problems have been shown in 
Figure (4,4) • It should be noted that the model STZDES-22 
is capable of designing the electrolyte hole profile for any 
desired spike profile. In the present work, the spike profile 
has been chosen such that it is similar in nature to the spike 
profile obtained by experimental tool, but lesser in height. 
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FIG. 4-3 DESIGNED TOOL FOR A STRAIGHT SIDED HOLE WITH A FLAT TOP SURFACE DURING ECD 
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FIG. 4.4 DESIGNED TOOL FOR A DESIRED SPIKE PROFILE IN STAGNATION ZONE DURING ECD 




(b) 

FIG. 4-5 Temperature, Conductivity and Void fraction 

variation along electrolyte flow path for Job No.427(EC 

(a) stagnation, front, transition 8c side zones 

(b) stray current zone- 
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CHAPTER V 

_C0NCLUSI0NS AND SCOPE FOR FUTURE WORK 

# 

5.1 CONaUSlONS 

The following conclusions have been drawn from the 
present work: 

- The model STZFET-22 is capable of predicting the 
nature of variation of current density, temperature,' 
electrolyte condbctivity etc, in all the five zones 
during Electrochemical drilling, i 

- Designed tool (cathode) shape is not the comple- 
mentary shape of desired workpiece shape reduced by 
some machining gap dimensions. 

There is a good agreement between the designed tool 
and the tool used during experimentation. Thus 
correction method is capable of designing tools 
reasonably accurately for Electrochemi.cal drilling. 

- The accuracy of correction method depends upon the 
accuracy of anode shape prediction model. Any 
discrepancy is anode shape prediction model will be 
reflected in the designed tool,. 

5.2 ■ SCOPE FOR FUTURE WORK 

The work embodied in this thesis can be extended in 


the following directions: 
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( i) An attempt should be made to improve the accuracy of 
the tool design model by improving the anode shape 
prediction model in the following ways: 

(a) More accurate electrolyte flow velocity distri- 
bution should be obtained by solving Navier- 
Stoke's equations, 

(b ) In the present analysis, effect of void fraction 
on conductivity in stagnation zone and stray ' 
current zone has been neglected. It would be of 
interest to know- to what degree this assumption caused 
the predictions to deviate from measured results. 

So effect of void fraction on conductivity must be 
analysed, 

( ii) Extensive data should be collected about the machining 
efficiency, electrochemical equivalent of different 
metals and alloys under different situations, so that 
accurate prediction of anode profiles can be made, 

( iii) Careful experiments should be conducted to find out 
the spike profile in stagnation zone. The final 
comparison between the experimental and calculated 
results must be left to show whether the assumptions 
made in the analytical model leaves the result suffi- 
ciently legitimate for the purpose in hand, 

( iv) For better accuracy, models based on two and three 
dimensional analysis of ECD using isoparametric and 
higher order elements should be evolved. 
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(v) For the practical applications (for example, 
elliptical, parabolic shape drilling, etCo) 

the tool design models based on three dimensional 
analysis should be evolved, 

(vi) This tool design procedure can be extended to 
external shaping as well^i. 


T 
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APPENDIX - I 




COMPUTATIONAL DETAILS 

This appendix describes the development of generalised 
computer programs for the computation of anode profiles obtained 
during ECD using a bare tool ( STZFET-22) and for computation 
of designed tools during ECD ( STZDES-22) In these models, 
three noded triangular elements have been used for discre- 
tizing the solution domaine The potential distribution has 
been obtained by solving the Laplace equation. These programs 
are capable of automatically generating the global stiffness 
matrix, applying the prescribed boundary conditions and storing 
the stiffness matrix in the banded form. The resultant set 
of equations have been solved by Gauss Elimination Technique, 

I 

The computational schemes used in anode shape prediction model 
STZFET-22 and tool design model STZDES-22 are shown in 
Figures A.1 and A. 2 respectivelyo 

BRIEF DESCRIPTION OF SUBROUTINES USED 

( I) SUBROUTINE INDATA: This subroutine reads and prints 

the input data for the problerao 

(II) SUBROUTINE XYCORD: It generates the node numbers, 

X and Y— coordinates for each node, connectivity matrix 
for the finite element mesh, the various arrays and 
assigns the boundary conditions for the nodes on the 


tool and workpiece^ 
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( III) SUBROUTINE BANDW: It determines the band width of the 
global stiffness, matrix, 

( IV) SUBROUTINE STIMA.T: This subroutine generates the 
coefficients of elemental stiffness matrices and 
stores the assembled global stiffness matrix in 
banded f orm, 

(V) SUBROUTINE ABCOND: This subroutine applies prescribed 
boundary conditions for the nodes on the tool and 

the workpiece, 

(VI) SUBROUTINE BANSOL: This solves the set of linear 

simultaneous equations by Gauss elimination method 
and gives the potential distribution, , 

(VII) SUBROUTINE VELSTI; This subroutine generates the 
elemental stiffness matrices and assemble them is global 
stiffness matrix. It also generates the velocity 
boundary matrix, 

(VIII) SUBROUTINE ABCVEL: This subroutine applies the velocity 
boundary conditions. 

(IX) SUBROUTINE VELSOL: This subroutine solves the set of 
linear simultaneous equations by Gauss’s elimination 
method, 

(X) SUBROUTINE VEL CAL: This subroutine calculates the 
velocities at various nodes on the workpiece, 

(XI) SUBROUTINE TEMCON: This subroutine calculates the 
temperature and current density at various nodes on 

the workpiece. 



(XII) SUBROUTINE COmOD: This subroutine calculates 

the electrolyte pressure and void fraction within 
the IB3. It also modifies the electrolyte conductivity 
by taking into consideration the void fraction* 

(XITI) SUBROUTINE FEDRAT: This subroutine determines the feed 
rate at various nodes on the workpiece* 

(XIV) SUBROUTINE NiRR: . This subroutine calculates the MRR 
and IBG at different nodes on the workpiece, 

(XV) SUBROUTINE MODY: This subroutine generates new elements 
at the end of each computational cycle and calculates 
the Y-coordinates of each node* 

(XVI) SUBROUTINE MODX: This subroutine generates the X-coor- 
dinates of each node at the end of each computational 
cycle, 

(XVII) SUBROUTINE PRTOUT: . This subroutine prints out the 
reguired results in a specific format, 

(XVIII) SUBROUTINE INITIL: This subroutine initializes the 

global stiffness matrices and the potential distribution 
matrices in each cycle of computation, 

SUBROUTINE DESIGN: This subroutine calculates the 
error and correction and develops the new tool shape,. 


(XIX) 




(a) 
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( START ) 



READ AKID PRINT INPUl' DATA (iNDATA) 


COMPUTE NO. OF COMPUTATIONAL CYCLES 


GENERATE X-COOR. and Y-COOR. FOR FEM MESH (XYCORD) 
^ — i ^ " 

DETERMINE BAND WIDTH (BANDW) 

- - -- -- 

GENERATE ELEMENTAL STIFFNESS MATRICES (STIMAT) 


I 



ASSEMBLE & STORE IN BA^DED FORM 

1 






APPLY BOUNDARY CONDITIONS (ABCOND) 

r— ' 



SOLVE SET OF EQUATIONS TO DETERMINE VOLTAGES (BANSOL) 





COMPUTE VELOCITIES AT VARIOUS NODES (VELCAL) 



1 





COMPUTE J,T AND K FOR Aa NODES (TEMCON) 








COMPUTE FEED RATE AT ALL NODES ON WORK (FEDRAT) 



V 




COMPUTE MRR & lEG AT ALL MODES (MRR) 


^ ^ ■ - 
PRINT RESULTS (PRTOUT) 

REASSIGN Y COORDINATES (MODY) 
REASSIGN X COORDINATES (MODX) 


No 



Yes 


STOP 


END 




FIGURE A.1 


(b) 

FLOW CHART FOR MODEL STZFET-22 

(a) FOR FRONT, TRANSITION, SIDE & STRAY CURRENT ZONE 

analysis a^tatvctc 

( b) FOR STAGNATION ZONE ANALYSIS 
















— 1 

SOLVE SET OF 

veL. potenti; 

1 f 

EQUATIONS TO GET 
\L (VELSOL) 





compute velocities at different 
nodes (VEL CAL) 





90 



(a) 
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(b) 


FIGURE A„2 ; FLOW CHART FOR MODEL STZDES 22 

(a) FOR FRONT, TRANSITION, SIDE & STRAY CURRENT 
ZONE AN/LYSIS 

(b) FOR STAGNATION ZONE ANALYSIS 
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APPENDIX II 
TABLE 1 


Identification Number and Hardness of Work Materials [3 ] 


Work Material 

Identification 

number 

Average 
hardness 
( BHN) 

Low alloy steel 

casti ngs 

200-299 

239 

Low alloy steel 

casti ngs 

400-499 

187 

Low alloy steel 

f orgings 

600-699 

21 8 

Low alloy steel 

f orgings 

700-799 

220 


Work Material Composition 

Alloying 

element 

Low alloy steel castings 

200-299 400-499 

C 

0,094 

0 o 234 

S 

0,018 

0.01 1 

P 

0.020 

0.014 

Si 

0,317 

0,500 

Mn 

0,400 

1 .440 

Ni 

- 

0.095 

Cr 

0,760 

0,260 

Mo 

0.920 

0.073 

Cu 

- 

0,065 

A1 ^ ^ 

- 

0 , i 055 

Va 

Very small 

- 

other impurities 

- 

5 to 6 


93 


TABLE 2 

Experimental Conditions used during ECM tests [ 3] 


Work material 


Tool material 

Tool diameter 

Mode of electrolyte 
supply 

Applied voltage 

Electrolyte flow 
ra te 

Electrolyte 
conductivi ty 


M.S. low alloy steel casting and low 
alloy steel f orgings 
Brass 

7 mm to 1 8 mm 
Radial outward flow 

t 

1 0-60 V 

6 —5 3 

6.7x1 0 ■■ to 15x1 0 m /sec 

-1 -1 

0o003 to 0„01 1 Ohm mm 


Machining time 


600 to 3600 sec 



TABLE 3 
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Experimental Conditions used 

durinq ECD of Cast 

Low Alloy Steel | 3 ] 


f = 0,000571 mm/s; Y. = 

1 ,0 mm 


Workpiece 

No. 

h 

( Ohm”* mm“^ ) 

E 

V 

( volt) 

( ampere) 

% 

Tool-1 

o 

9 

II 

(N 

rtc=1 .246 



413 

0,. 00725 

17,75 

13.1 0 

90.15 

408 

0„00725 

15.31 

16.24 

91 .70 

409 

0.00725 

21 .84 

14.14 

95.00 

406 

0.00725 

12.94 

13.28 

86.90 

Tool-3 

r^ = 4,15 




407 

0.0063 

12.93 

15.00 

93.15 

421 

0.0070 

14.89 

16,40 

96.20 

424 

0,i0060 

18,50 

19.86 

92.00 

425 

0,0066 

20;i49 

18.80 

91 .86 

Tool-4 

= 4,505 

rtc=2,13 



429 

0,0060 

11 ,57 

1 8.53 

94.00 

401 

0.00547 

1 2.975 

22.15 

91 .82 

41 8 

0,00532 

15.95 

18.40 

92.50 

427 

0.00532 

17,i62 

— 

— 

402 

0.00600 

15,87 

1 8.40 

93.70 
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TABLE 4 


Voltage Variatio n during Electro-chemical Drilling (bare 


type oF tool) [3] 


Job No, 


Time in 
min. 

421 

425 

407 

409 

424 

408 

427 

0 

20.„5 

21 .25 

19,;5 

20.75 

21 ,0 

16,0 

24.0 

1 

20,0 

21 .0 

18,5 

20.25 

20.5 

15.5 

23,0 

2 

19,0 

20.75 

17.5 

20,0 

19.0 

1 8.0 

22.0 

3 

18,0 

20.25 

17.0 

20.0 

1 8.0 

17,5 

21 ,0 

4 

18,0 

20.0 

15.5 

19,8 

18.25 

17.5 

20,5 

5 

16,5 

20.0 

14.5 

19.75 

17.75 

16.5 

19,0 

6 

15,5 

19.8 

14,5 

19.8 

1 8.0 

1 6,5 

1 8.75 

7 

15.5 

19.75 

14.0 

19,8 

1 8,25 

16,5 

1 8.5 

8 

15,0 

19,8 

14.0 

19.6 

18,5 

16,25 

18.5 

9 

15,0 

19,8 

13,75 

19,25 

1 8.25 

16.0 

17,0 

1 0 

15,0 

19.6 

13.0 

19,0 

18,0 

15,75 

17.0 

1 1 

14.5 

19,25 

12.5 

19.0 

18.0 

15,40 

16,5 

12 

14,0 

19.0 

12,5 

19,0 

17.75 

15.0 

16.5 

13 

14.5 

19.0 

12.0 

1 8;;8 

17.5 

15.0 

16,5 

14 

14.5 

19,0 

11 ,75 

18.5 

17,5 

15.0 

16,5 

15 

14.5 

18,8 

12.0 

1 a ,5 

17.5 

14.75 

16.625 

16 

14.0 

18,5 

1 1 ,5 

1 8„5 

17,0 

14,75 

16,5 

17 

14,0 

18,5 

11.5 

18.5 

16.875 

15,0 

16,5 

18 

13.5 

18,5 

11 ,0 

1 8,25 

15.875 

15,5 

16 i 5 

19 

13.5 

1 8, ,5 

11,0 

1 8.5 

16,5 

15.5 

16,5 

20 

13.5 

18.5 

11 .0 

18.5 

16.5 

16.0 

16,5 

21 

13.5 

18, ,5 

11,0 

1 8;,5 

16.0 

17,0 

16.5 

22 

13.5 

18.5 

1 0.75 

18 i .5 

16.5 

17.0 

16,5 

23 

13,5 

1 8,5 

10,75 

1 8..5 

15.5 

16.0 

16.5 

24 

13.5 

18,5 

11 ,25 

18.5 

15.5 

16,0 

16,5 

25 

13,25 

1 8.5 

11 ,0 

1 8,25 

1 5 .5 

13,0 

16,0 

26 

13,0 

18.5 

11 .0 

1 8,25 

15.5 

13,0 

16,0 

27 

1 3.5 

1 8,5 

11,0 

1 8.25 

1 5 .5 

13,0 

15.75 

28 

13.25 

18,25 

1 1 ,0 

1 8.0 

15.5 

13.0 

1 5 .75 

29 

13.0 

1 8.25 

11.0 

1 8.0 

1 5 .5 

13.0 

15,5 

30 

13.0 

18,0 

11 ,0 

18.,0 

15.5 

13.0 

15,5 
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TAELE 5 



Job Noo 

425 

427 

421 

409 

408 

406 

Depth 

X,mm 



Overcut, 

£nm> 



0,0 

. 2.26 

2,32 

1.71 

2,45 

2.634 

1 .902 

0,5 

1 .76 

1 ,73 

1 .54 

2.01 

2,46 

1 .52 - 

1.0 

1 .695 

1.65 

1.56 

1,84 

2.214 

1 .75 

1.5 

1 .61 

1 ,63 

1.4 

1.7 

2,14 

1 ,33 

2.0 

1 .6 

1 .44 

1 .45 

1 .67 

2,02 

1.36 

2,5 

1 .51 

1 .51 

1 .31 

1 .57 

1 .964 

1,21 

3.0 

1 .52 

1 ,34 

1 .36 

1 .52 

1 .85 

1 ,25 

3.5 

- 

1 .42 

. 1.13 

- 

1.76 

1.11 

4.0 

1.43 

1 .12 

1 ,25 

1 .42 

1 .67 

1,154 

4,5 

1 ,33 

1,31 

1 ,02 

1 ,38 

1 .57 

1,03 

5.0 

1.34 

1,02 

- 

1.34 

1 ,48 

1,05 

5,5 

1.23 

1.18 

0,.91 


1.34 

0,924 

6.0 

1 .11 

0.91 

0.92 

1.26 

1.28 

0.96 

6.5 

1 .1 

1 ,03 

0.71 

1 .26 

1.17 

0,84 

7.0 

0.91 

0,61 

0.715 

1 .16 

1 .134 

0,87 

7,5 

0,84 

0,61 

- 

1 .11 

0.95 

0.74 

8,0 

•> 

- 

0.60 

0^82 

0.78 

0,80 


